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Abstract - Evidence from NMR of a two-component spin system in cuprate high-T c super- 
conductors is shown to be paralleled by similar evidence from the electronic entropy so that a 
two-component quasiparticle fluid is implicated. We propose that this two-component scenario is 
restricted to the optimal and underdoped regimes and arises from the upper and lower branches 
of the reconstructed energy-momentum dispersion proposed by Yang, Rice and Zhang (YRZ) to 
describe the pseudogap. We calculate the spin susceptibility within the YRZ formalism and show 
that the doping and temperature dependence reproduces the experimental data for the cuprates. 



From the electronic entropy we present evidence for 
a two-component electron fluid in underdoped high-T c 
cuprates matching similar evidence from NMR p}[3]. We 
then show that this two-component behavior probably 
arises from Fermi surface reconstruction and the associ- 
ated band-splitting that occurs due to pseudogap corre- 
lations. We illustrate this using the model proposed by 
Yang. Rice and Zhang (YRZ) [3]. Implicit in this inter- 
pretation is the prediction that single-component behavior 
is recovered in the overdoped region where the pseudogap 
closes around p » 0.19 holes/Cu. 

Hole-doped HTS, in their overdoped state, possess a 
large Fermi surface [MZ] enclosing (1 + p) carriers where 
p is the doped hole concentration residing largely on in- 
plane oxygen orbitals and the "1" arises from the unpaired 
electrons residing on the Cu sites. This naively suggests 
a two-component electron system. However, Zhang and 
Rice [5] showed that the doped holes form a local sin- 
glet state with the hole on the Cu 3d x 2 _ y i orbital, the 
so-called Zhang-Rice singlet. It was suggested that the 
singlet binding energy was so large that the oxygen 2p 
holes do not contribute to the spin susceptibility thus re- 
sulting in a single-component spin scenario. This situ- 
ation has been considered long-established starting from 
Takigawa et al. [9], who showed that the 17 and 63 Cu 
Knight shifts, 17 K(T) and 63 K(T) respectively, for un- 



derdoped YBa2Cu307_,5 exhibit an identical temperature 
dependence. Subsequent 89 Y NMR studies confirmed the 
same T-dependence in 89 K(T) [TU] adding further weight 
to the single-component scenario. 

However, this single-component picture was recently 
questioned by Haasc et al 1 3_ . In order to eliminate the 
Meissner term in the spin shift they constructed a term 
Gj_ by subtracting the Knight shift for the apical oxygen 
with field perpendicular to the c-axis from that for the 
planar Cu shift, and similarly G\\ for field parallel to the 
c-axis. Thus, 



G±(T) 

G||(T) 



63. _ 



K{T) - l7 ' A > x K(T), 



17,P,| 



\K(T) 



_17,A,|| 



K(T), 



(1) 



where the superscript A refers to the apical oxygen while 
the superscript P refers to planar oxygen. 

These authors showed that G±(T) displayed a T- 
independent Pauli-likc metallic susceptibility, while G||(t) 
displayed a strongly T-depcndcnt susceptibility consistent 
with a gapped normal-state spectrum. By constructing a 
two-component ansatz: 



G±(T) =cnxi +C12X2, 
G\\(T) = C21X1 + C22X2, 



(2) 
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Fig. 1: (Color online) (a) Reproduced from ref. [9]: the Knight shift, K(P), for YBa 2 Cu 4 08 at P = 20, 42 and 63 kbar plotted 
vs K(l bar) where T is the implicit variable. Arrows indicate T c . Above the range of SC fluctuations K(P) is linear in K(l 
bar), (b) An analogous plot of S(T,p)/T for YBa 2 Cu3 06+z plotted vs S(T,p = 0.133)/T for x = 0.76, 0.80, 0.87, 0.92 and 
0.97 (open circles). The corresponding p values are listed. Also shown is RwXs(T,p) vs RwXs(T,po — 0.135), expressed in 
entropy units (solid curves), (c) S(T,p)/T for Bi2Sr2CaCu2 0s+5 plotted vs S(T,po — 0.129)/T for the listed p-values. 



where \i an d Xi are the two putative uniform susceptibil- 
ities, Haase et al. [Tj showed that Xi(^) has a characteris- 
tic pseudogap-like T-dependence, falling steadily towards 
zero with reducing T, while X2(T) has a constant Pauli-like 
behavior and only falls to zero below T c with the opening 
of the superconducting (SC) energy gap. 

In introducing xi an< i X2 Haase et al. took a 
lead from Johnston |11) who, from static bulk measure- 
ments, inferred a two-component susceptibility of the form 
X(p, T) = Xi{Pi T) + X2(p) where %2 is a function of dop- 
ing only. However, the features he sought to explain arise 
naturally from the presence of the van Hove singularity 
(vHs) in the overdoped region which ensures an overall 
increasing density of states across most of the phase dia- 
gram (as reflected in the X2(p) term). A single-component 
susceptibility fully accounts for the complete evolution 
of S9 K S across the undcrdoped and overdoped regions in 
Yo.8Gao.2Ba 2 Cu307-5 [12] provided the full ARPES de- 
rived dispersion is utilized, including the vHs. It requires a 
differential measurement, as in Eq.(l), to expose the true 
underlying two-component behavior. 

To make the problem more specific we consider the most 
recent study from the Haase group [3], which also sup- 
ports the two-component picture. In this they carried 
out high-pressure diamond-anvil NMR measurements on 
YBa2Cu408 to 63 kbar. Pressure, P, acts dominantly, 
though not exclusively, to increase the doping by trans- 
ferring holes from the CU2O2 chains to the CuC>2 planes. 
Thus one typically observes a pressure-induced decrease 
in thermoelectric power [13] consistent with an increase 
in hole concentration, p, resulting in a pressure-induced 
increase in T c on the undcrdoped side and a decrease in 
T c on the overdoped side [M]. (That the effect of pres- 
sure is not merely to increase the hole concentration is 
evident from the fact that the maximum T c is increased 
from approx 93 K at ambient pressure to approximately 



107 K at a pressure of about 7 GPa [E]). Consistent with 
this, Meissner et al. [3] observed 17 K(T) to progress from 
a strongly T-dcpcndcnt behavior under ambient pressure, 
typical of a pseudogapped underdoped cuprate, towards a 
nearly T-independent behavior at a pressure of 63 kbar, 
more typical of an optimally-doped cuprate. By plotting 
17 K(T, P) versus 17 K (T, P = 1 bar) for P = 20, 42 and 63 
kbar with T as the implicit variable they obtained a lin- 
ear relation that showed a characteristic progression with 
increasing P. This plot is reproduced in Fig. QJa). 

The authors draw attention to the fact that well above 
T c there is a linear region that for the lowest doping ex- 
tends down almost to T c . They model this linear behavior 
within a two-component scenario. 

Firstly we wish to point out that identical behavior is 
seen in the electronic entropy, S(T), thus indicating that 
it is not just the spin system that exhibits two components 
but the total quasiparticle ensemble. 

For a nearly free-electron system the spin susceptibility, 
Xs and S/T are related via the Wilson ratio, Rw'- 



RwXs{T) = S(T)/T. 



(3) 



Accordingly, the open circles in Fig. [IJb) and (c) show 
S(T,p)/T versus S(T,p )/T (where T is the implicit vari- 
able) for YBa2Cu306+x (b) and for Bi2Sr2CaCu20s+<5 (c). 
We choose po = 0.133 for the former and po = 0.129 
for the latter, which are close to the zero-pressure doping 
state of YBa2Cu408, po ~ 0.13. The same generic behav- 
ior also occurs in plots of S(T,p)/T versus S(T,po)/T for 
Yo.8Cao.2Ba 2 Cu306+:r (not shown). The correspondence 
between K(P) and S(p)/T is remarkable. 

Assuming the major effect of pressure is to in- 
crease doping and letting a be the thermopower, then 
(da/dP) T=2 9o = -0.075^V K^kbar- 1 [13]. Using the 
thermopower correlation with doping (da/dp)T=290 = 
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-134^V K^hole" 1 US one finds (dp/dP) = 5.6 x l(r 4 
holes/kbar. Thus the doping states of YBa 2 Cu40s at 20, 
42 and 63 kbar are 0.141, 0.154 and 0.165 holes/Cu. Vi- 
sually the red, green and blue data sets in Fig. 0Ja) cor- 
respond closely to the red, green and blue data sets in 
Fig. HJb) an d so should be at roughly the same doping 
states. This is indeed the case where in Fig. [TJb) the dop- 
ings are seen to be 0.140, 0.148 and 0.163 holes/Cu. 

Moreover, we show that the correspondence between 
Xs and S/T is in excellent quantitative agreement with 
Eq. [3l The solid curves in Fig. [ljb) show values of 
Xs{T,p) versus Xs(T,po = 0.135) expressed in entropy 
units using the Wilson ratio for nearly-free electrons R^ = 
(7r 2 /3)(fcs//Us) 2 . We use the bulk susceptibility Xs data 
of Cooper and Loram [16]. One can see that the cor- 
respondence between Xs and S/T is not just qualitative 
but quantitative. These plots reveal precisely the same 
magnitudes as S/T and the same breakaway from linear 
behavior sets in well above T c due to strong SC fluctua- 
tions il7 . (Note that below T c the diamagnetism leads 
to different behavior from S/T). We thus conclude that 
there is in fact a two-component quasiparticle system, not 
just a two-component spin system. 

We now turn to our central thesis that this generic be- 
havior, seen in both the spin susceptibility and the en- 
tropy, arises from band splitting that occurs when the 
Fermi surface reconstructs due to competing orders such 
as a CDW, SDW or short-range AF correlations. We 
illustrate this within the YRZ scenario where the elec- 
tron self-energy term E^(k)/(uj + ££) reconstructs the 
E(k) dispersion into upper and lower branches which 
yield our two-component quasiparticle ensemble. Here 
££ = — 2t(p)(cosk x + cos k y ) is the nearest- neighbor term 
in the tight-binding E(k) dispersion. The coherent part 
of the electron Green's function is given by: 



G(\a,w,p) = gt(p) 



6c- 






(4) 



where £& = — 2t(p)(cos k x + cos k y ) — 4t'(p) cos k x cos k y — 
2t"(p)(cos2k x + cos2k y ) — fj, p (p) is the tight-binding 
dispersion to third- nearest neighbors, and E g (k) = 
^Eg(p)(cosk x — cos k y ) is the pseudogap with doping de- 
pendence Eg(p) = 3£o(0.2 — p), while for p > 0.2 we 
have Eg(p) = 0. (This means the pseudogap closes at 
p = 0.2 however we have extensively shown this to oc- 
cur at slightly lower doping p = 0.19 [18j . Here we re- 
tain the value 0.2 to remain consistent with YRZ). The 
chemical potential fi p {p) is chosen according to the Lut- 
tingcr sum rule. The doping dependent coefficients are 
given by t(p) = g t (p)t + (3/8)g s (p)Jx, t'(p) = g t {p)t' 
and t"{p) = g t (p)t'o, where g t (p) = 2p/(l + p) and 
g s (p) = 4/(1 + p) 2 are the Gutzwiller factors. The bare 
parameters t' /to = —0.3, t" /to = 0.2, J/to = 1/3 and 
X = 0.338 arc the same as used previously [3]. 



Equation 2] can be re- written as 

9t{p)W£{p) 
w-E£(p) ■ 



G(k, U ,p)=Yl 



(5) 



a=± 



where the energy-momentum dispersion is reconstructed 
by the pseudogap into upper and lower branches 



1 



k 2 ^- «;-ni 2 



e£ = ^-&)-. 

which are weighted by 



-#'l ^(k), (6) 



W, 



1± 



(a + e£)/2 



v/[(a+^)/2] 2 + ^ 2 (k) 



(7) 



The spectral function is given by 

A(k, U ,p) - Y, 9t(j?WZ5{u - E£), (8) 



a=± 



from which the density of states can be calculated 

JV(w) = X;^(k»w)- (9) 



Finally, the spin susceptibility is given by 



X S = 2 M | J Y-^)jV(w)dw. 



(10) 



Figurc[2ja) shows the partial densities of states (PDOS) 
for p = 0.14 and 0.16 calculated from the upper and lower 
branches of the reconstructed dispersion given by Eq. [51 
The PDOS of the lower branch is quasi-linear across the 
Fermi level (w = 0) resulting in a roughly T-independent 
contribution to the susceptibility at low doping, shown in 
Fig. [3|b) by the solid curves. For dopings below about 
0.18, the PDOS of the upper branch lies above the Fermi 
level resulting in a gapped spectrum. With decreasing 
doping the upper PDOS is pushed further from u = 0, 
producing a T-dcpcndcnt contribution to x that is char- 
acteristic of the pseudogap state, shown in Fig. [5Jb) by 
the dashed curves. 

The two sets of susceptibilities shown by the dashed 
and solid curves in Fig. [5Jb) are to be directly compared 
with the functions Xi(T) and X2(T), respectively, reported 
by Haase et al., [TJ[2]. They reveal the same qualitative 
behavior but there are points of difference, primarily in 
their relative magnitudes where, at high temperature, the 
lower-branch susceptibility is more than three times the 
magnitude of the upper-branch susceptibility. By con- 
trast, Haase et al. find xi exceeds xi by around 50% 
at 300 K. We return to this discrepancy below. 

In order to compare directly with the pressure- and 
doping-dependent data shown in Figure Q] we plot in 
Fig. E^c) the sum of the susceptibilities (from the up- 
per and lower branches) as a function of the susceptibil- 
ity sum for p = 0.13, with T as the implicit variable. 



p-3 



Storey and Tallon 



o 0.2 



w 
o 

Q 





11 ' 




p 

0.16 
0.14 


J 




1 
I 

5 




' 


v^ V 4-"' 






0.0 0.2 0.4 

CO/t 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 

KJ/L 




0.07 0.08 0.09 0.10 

X(T, p=0.13) 



Fig. 2: (Color online) (a) Partial density of states calculated from the lower (solid) and upper (dashed) branches of the YRZ 
reconstructed dispersion, (b) Susceptibility calculated from the lower (solid) and upper (dashed) partial density of states, (c) 
Total Xs(T,p) plotted versus Xa(T,po = 0.13) where T is the implicit variable (units of \ are the same as in (b)). 



The qualitative features seen in the normal-state Knight 
shift (Fig. [T][a)) and electronic entropy (Figs. QJb) & ( c )) 
are reproduced in detail. In particular, the high-T down- 
turn seen for p = 0.17 reflects the proximity of the vHs. 
The same downturn is seen in the E^S^CaC^Os+s data 
(Fig. [BJe)) where the vHs indeed lies nearby [19] in the 
moderately overdoped region, whereas it is not evident in 
the data for YI^OusC^-^ (Fig. [ljb)) where the vHs is 
known to lie in the more deeply overdoped region. 

We are however left with two remaining questions: (i) 
why does G±(T) reflect more the lower branch suscepti- 
bility while G||(T) reflects more the upper branch suscep- 
tibility? And (ii) there is the question, mentioned above, 
of the relative magnitudes of the two susceptibilities. We 
suggest these have a common origin, as follows: 

The apical oxygen is coupled to the planar copper and 
oxygen orbitals via the Cu 4s orbital [5D] . This introduces 
matrix elements that weight the k-space sums, minimizing 
the contribution along the zone diagonals and maximizing 
contributions at the (tt,0) zone boundaries. This clearly 
will diminish the contribution to the susceptibility from 
the lower branch, effectively reducing the magnitude of the 
coefficients C\2 and cii- As a consequence the relative con- 
tributions of xi and Xi to G±(T) and G\\(T) differ, with 
G±(T) dominated more by xi- The c-axis hopping matrix 
is ij_(k) = txMk w here /j,^ = 0.5wk (cosk x — cosk y ). We 
find this does indeed bring the magnitudes of xi an d Xi 
closer together for lower doping but not so much at higher 
doping close to where the pseudogap closes. We await a 
more rigorous treatment of the precise interaction of the 
apical nucleus with the two global spin susceptibilities. 

In conclusion, we show that the entropy term S(T)/T 
displays the same doping evolution as the pressure- 
dependent Knight shift, thus indicating that the two- 
component electronic behavior resides in the quasiparticle 
spectrum and not just in the spin spectrum. We then show 
that the essential features of the two-component system 
are likely to arise from band splitting due to zone-folding 



effects as described e.g. by the Yang- Rice- Zhang model for 
the pseudogap. The pseudogap-like susceptibility xi in- 
ferred by Haase et al. arises from the upper branch and the 
Pauli-likc xi arises from the lower branch. If correct then 
it follows that single-component electronic behavior will 
be recovered when the pseudogap closes. Measurements 
of G± and G\\ will then help settle the still contentious 
issue as to whether the pseudogap closes abruptly, both 
above the SC dome and below it at a putative quantum 
critical point at p « 0.19. 
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